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ABSTRACT 

We discuss the peculiar nature of the nucleus of M94 (NGC 4736) in the context of new measurements 
of the broad Ha emission from HST-STIS observations. We show that this component is unambiguously 
associated with the high-resolution X-ray, radio, and variable UV sources detected at the optical nucleus of this 
galaxy. These multi-wavelength observations suggest that NGC 4736 is one of the least luminous broad-line 
(type 1) LINERs, with Lboi = 2.5 x 10'"' erg s"L This LINER galaxy has also possibly the least luminous 
broad line region known (Lna = 2.2 x 10^^ erg s~'). We compare black hole mass estimates of this system to 
the recently measured ^ 7 x 10* Mq dynamical black hole mass measurement. The fundamental plane and 
M-a* relationship roughly agree with the measured black hole mass, while other accretion based estimates 
(the M-FWHM{Ha) relation, empirical correlation of BH mass with high-ionization mid IR emission lines, 
and the X-ray excess variance) provide much lower estimates (^ IO^Mq). An energy budget test shows that the 
AGN in this system may be deficient in ionizing radiation relative to the observed emission-line activity. This 
deficiency may result from source variability or the superposition of multiple sources including supernovae. 



Subject headings: galaxies: active - galaxies: nuclei 
4736) 



1. INTRODUCTION: LINERS AND M94 

Most high mass galaxies are known to host massive black 
holes (BH), some passively lurking in their centers while oth- 
ers are actively accreting surrounding material (e.g., Giiltekin 
et al. 2009a and references therein). The mechanism that 
causes a BH's activity to turn on and off is still largely un- 
known. Understanding the structure of the active galactic nu- 
clei (AGN) at their lowest luminosities is crucial to determin- 
ing the physical and possibly evolutionary links between the 
most luminous galaxy centers and the passive ones. However, 
at low luminosities it is difficult to disentangle the various 
emission mechanisms that could be concurrently present in 
galaxy centers. As a consequence, the dominant power source 
of a large majority of actively line-emitting galaxies remains 
ambiguous (e.g.. Ho 2008 for a review). 

Diagnostic diagrams de .g.. Baldwin, Phillips, & Terlevich 
198 1\ 'Veilleux & O sterbrock ..1987: .Kauffmann et al. , ,2003i: 
Kewley et al. 2006) are relatively successful in separating 
out bona-fide accretion sources (Seyferts) from nuclei whose 
emission-line activity is mainly powered by young, hot stars 
(H II galaxies), based on emission line ratios. At least 
50% of the strong line emitters fall easily onto the H II lo- 
cus, however, only less than 10% are of the Seyfert type 
jHo, FiliDpenko & Sai-gent' '1997^ iKauffmann et al. I l2003t 
IConstantin & Vogeley 2006). A large fraction of the objects 
"in between" these two categories, that exhibit relatively low 
levels of ionization (i.e., low values of [O III]/H/3), maintain 
reasonably strong forbidden line activity (i.e., high values of 

' Based on observations made with the NASA/ESA Hubble Space Tele- 
scope, obtained from the data archive at the Space Telescope Science In- 
stitute. STScI is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. 

^ Some observations reported here were obtained at the MMT Obser- 
vatory, a joint facility of the Smithsonian Institution and the University of 
Arizona. 

' James Madison University, Harrisonburg, VA 22807 
University of Utah, Salt Lake City 841 12 
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[N II]/Ha), and are classified as low ionization nuclear emis- 
sion regions (LINERs). The reminder are usually called Tran- 
sition objects (Ts). Whether Ts and Ls are powered, at least 
partly, by accreting BHs, and thus could be called AGN, is a 
matter of continuous debate (Ho 2008). 

There are some typical emission characteristics that are 
considered to be particularly good indications that accretion 
onto a massive BH is an important, if not the dominant source 
of ionization in some Ls, and possibly Ts. The detection of 
broad Ha emission, regardless of its strength or luminosity, is 
generally considered to be the clue to AGN emission. Some 
LINERs (and maybe Ts as well) exhibit these features, how- 
ever, the majority of these systems show only narrow emis- 
sion, which could be generated by shocks, post-starbursts, or 
other processes unrelated with accretion. Observations out- 
side of the optical wavelengths often reveal AGN signatures in 
ambiguous and even starburst galaxy nuclei. X-rays are par- 
ticularly good tracers of accretion, however, they are not effi- 
cient in distinguishing AGN at Lx < 10'*^ erg s"', where con- 
tamination by X-ray binaries can be significant. X-rays are 
also unlikely to detect heavily absorbed AGN (i.e., Compton- 
thick; A^H > 1 -5 X 10^"* cm"^). On the other hand, mid-IR high 
ionization emission lines like [Ne V] A14.32 /im, 24.32/im 
(97.1 eV) appear to be a trustworthy indicator of AGN activity 
dWeedman et al. IIIOOllArmus et al. ll2006HAbel & Satvapall 



"2008^ due to the extreme conditions (i.e., very hard ioniz 
ing radiation) required to produce them; because these fea 
tures have considerably lower optic al depth, their dete ction 
can also reveal Compton thick AGN (ISturm et al. 1120021 e.g., 
NGC 1068). This technique has now been appUed to re- 
veal new and la rge numbers of optically unidentified AGN 
dSat y apal et"an r2008: Goulding & Alexander 2009), provid- 
ing thus sensitive improvements on previous AGN censuses. 
These studies, along with X-ray an d radio studies of nearby 
galax i es (e.g. iN agar. Falcke, & Wilign] l200l IZhanget al l 
l2009t iDesrochTs. Greene & Ho ..2009^ suggest that a major- 
ity of LINERs and a large fraction of transition galaxies might 
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in fact host accreting black holes. The presence of an accret- 
ing black hole does not however guarantee that the accretion 
power is the dominant source of ionization of those galaxy 
nuclei. A more recent assessment of the energ y budget of 
LINERs bv'Eracleous. Hwang. & F lohic I (1201 Obi) argues that 
in 85% of LINERs the AGN ionizing photons are not suffi- 
cient for producing the observed nebular emission, and thus 
other power sources are likely to dominate. 

Some new potentially powerful insights into the excitation 
mechanism of the low luminosity AGN (LLAGN), and in par- 
ticular the ambiguous sources, come from recent studies of 
large statistical samples of nearby galaxy nuclei, which reveal 
a potential // II — > 5 — > T -> LINER — > Passive Galaxies evo- 
lutio nary sequence in the process of BH gr o wth within galax- 
ies (iConstantin. Hovle. & Vogelev 1 120081 : IConstantin et al. I 
I2009h . This sequence traces trends in (1) increasing host 
halo mass, (2) increasing environmental density, (3) increas- 
ing central BH mass and host stellar mass, (4) decreasing BH 
accretion rate, (5) aging of the stellar population associated 
with their nuclei, and (6) decreasing in the amount of dust ob- 
scuration, which might translate into a decrease in the amount 
of material available for star-forming or accretion. In this pic- 
ture, Seyferts and Ts are transition phases between the initial 
onset of accretion, usually swamped by the star-forming gas 
and associated dust, which is seen optically as an H II sys- 
tem, and the final phase of accretion observed as LINERs of 
already massive BHs. While this idea is supported by vari- 
ous other independent observ ational studies of low luminosit y 
AGN and starburst galaxies (ISchaw inski et al. "2007. ' 20 1 Oh . 
along with state of the art hydrodynamical models for the life- 
cycles of the most luminous AGN (e.g., Di Matteo et al. 2005, 
Hopkins et al. 2006), it is very probable that not all sources 
fit into this scenario. It is also important that this evolutionary 
sequence idea is tested on samples that span a narrow distribu- 
tion in Hubble types; because most of the H lis and Ts are rel- 
atively late-type disk galaxies that, likely, never experienced 
a (recent) major merger, it is possible that the trigger of such 
a sequence is different from the merger that initiates a simi- 
lar life-cycle at high luminosities. Nevertheless, investigating 
the challenges certain objects bring to this idea is useful for 
identifying and quantifying the caveats associated with this 
sequence. These objects may also be at interesting stages in 
their galaxy evolution. 

NGC 4736 (or M94, UGC 7996) is a captivating example 
of an ambiguous galaxy nucleus, which poses challenges to 
the general understanding of AGN phenomena, including the 
above mentioned sequence. This object is one of the clos- 
est {d = 4.3 Mpc; Tully & Fisher 1988) nearly face-on spi- 
ral, with a SAab Hubble Type. Its proximity enables study 
of details that would be unobservable in more distant sys- 
tems. Its nucleus has a low-luminosity LINER spectrum, 
but has been also included in catalogs of Transition objects 
dFTlippenko & Sargent! 119851) . or Seyfert 2s (e.g., SpinelH 
et al. 2006). The AGN nature of this object has been 
constantly debated, an aging starburst being a compelling 
alternative jMaoz et al. |[l99l l200l lEracleous et adl2002[ 
iKording et al. 112005b . The galaxy presents a ring of H II re- 
gions at a radius of ^ 50", red arcs at ^ 15", a high surface 
brigh tness nuclear region, and high far - infrar ed bulge emis- 
sion dKinnev et al. 1 119931: iSmith et al~l 119941) . Its intricate 
structure of off-nuclear comp act source detections in X-ray 
(lEracleous et al. 2002), radio ( IKording et al. II2005I) and UV 
dMaoz et al. Ii200 5). that do not necessarily match with each 
other, certainly increase the ambiguity associated with the na- 



ture of the main nuclear ionization mechanism. A common 
implied scenario in all of these studies is that this system is 
probably in the final stages of a merger 

We reexamine here NGC 4736 in the context of additional 
evidence for its AGN nature, which is the detection of a 
broad Ha component in its nuclear spectrum, as observed 
by the Hubble Space Telescope {HST) with the Space Tele- 
scope Imaging Spectrograph (STIS). We also gather multi- 
wavelength data and show that source of the broad H-alpha 
emission line is coincident with a compact X-ray and radio 
source. These observations suggest that NGC 4736 hosts 
a broad line region of significantly low luminosity, which 
makes this object one of the least luminous LINERs with 
strong evidence for BH accretion. Interestingly, BH mass in- 
dicators calibrated on rapidly accreting Seyfert galaxies give 
highly discrepant mass estimates in this more quiescent sys- 
tem. The ^ 1O^M0 value given by the M-a* is two orders 
of magnitude higher than the values obtained via estimators 
based on the observed emission (X-ray variability, scaling re- 
lations, mid-IR emission) which, although physically inde- 
pendent of each other, give a consistent result of ^ 1O^M0. 
We are thus facing the following conundrum: either a) the 
standard AGN BH mass indicators do not necessarily apply 
to sources emitting in this low luminosity regime as, proba- 
bly, the emission mechanism is fundamentally different from 
that associated with higher luminosity AGN, or b) the emis- 
sion signatures do not trace accretion onto the central BH. We 
discuss possible resolutions of this discrepancy in §5, where 
we propose some rather exotic scenarios. 

2. DATA COMPILATION & ANALYSIS 

NGC 4736 has been quite extensively observed across the 
whole electromagnetic spectrum. We present in this section 
the multiwavelength observations of this galaxy nucleus, in 
connection to new measurements of the broad Ha emission 
detected with HST -STIS. 

2.1. The broad Ha emission observed with HST-STIS 

High resolution optical spectra of the NGC 4736 nucleus 
were obtained with HST-STIS on July 20020 The observa- 
tions were obtained with the 52" x 0"1 aperture oriented at 
PA= 49.65°, with the slit centered along the major axis of the 
starlight distribution; two cosmic ray split exposures were ob- 
tained, one being slightly shifted in the slit direction. The 
total combined exposure time is ^4000s. The G750M grating 
was set at 6581 A, with a scale of 0"05/pixel, with no bin- 
ning. We reduced the spectra using IRAIQ and the STIS re- 
duction pipeline maintained b y the Space Telescope Science 
Institute (D ressel et al. I20()7h . This reduction included image 
combination and cosmic-ray rejection, flux calibration and 
correction of the wavelength to the heliocentric frame. To 
measure the nuclear nebular line-emission properties we ex- 
tracted the 1 -dimensional aperture spectrum five pixels wide 
(0."25) centered on the continuum peak. The extracted spec- 
trum thus consists of the central emission convolved with the 
STIS spatial point-spread function (PSF) and sampled over a 
rectangular aperture of 0"25 x 0"1. The measurements span 
6295 - 6867A with a resolution of 0.87A (cri„st = 17 km s"'). 

Figure[T]shows the resulting spectrum together with the best 
matching continuum model, spectral fits of the pure emission- 

^ Data are publicly available (Prop ID 8591) but have not been published. 
* IRAF is distributed by NOAO, which is operated by AURA Inc., under 
contract with the National Science Foundation 



The M 94 conundrum 



3 



TABLE 1 

Broad Ha Emission Measurements 
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Fig. 1.— Spectral fits to the HST-STIS spectrum of NGC 4736, based on 
a 5 pixel extraction aperture, which corresponds to a total emission area of 
0.025 arcsec- . Upper panel: The thin continuous line shows the observed 
data, while the green line shows the best fit stellar population model obtained 
by clipping out the hashed emission-line region. Middle panel: The contin- 
uum subtracted spectrum, i.e., the emission-line component, together with 
the corresponding spectral fit, that have been shifted by an additive constant 
to match the initial continuum level, for illustration purposes. The thin con- 
tinuous lines show the individual Gaussian components ( 1 per narrow line), 
the thick grey continuous line represents the broad Ha feature, and the red 
line is the final fit to the observed spectrum. Lower panel: The residuals after 
the subtraction of the model fit. 

line component, and the associated residuals. The continuum 
model of the underlying stellar population is obtained via a 
minimization of a nonnegative least-squares fit between the 
observed spectrum and a s um of discrete star bu rsts of differ- 
ent ages, adopted from the iBruzual & Chariot I (^03) stellar 
population synthesis templates, together with dust attenuation 
modeled as an additional free parameter. The continuum fit- 
ting is performed using an adaptation to our data of Christy 
Tremonti's code (Tremonti et al. 2004). The pure emission- 
line spectrum, obtained by subtracting the modeled contin- 
uum from the observed spectrum, is fit by a combination of 
linear continuum and Gaussian components. In addition to 
narrow emission, the Hq!-i-[N II] feature shows clear evidence 
for a broader feature. A flux ratio of 1:3 was assu med for the 
[N II] doublet, as dictated by the branching ratio (OsterbrockJ 
Il989); the [O I] feature and the [N II] and [S II] doublets were 
assumed to share a common velocity (red)shifts and widths. 
The best-fitting Gaussian parameters were derived via an in- 
teractive minimization, using SPECFIT (iKriss II 19941) . Be- 
cause of the generally low signal to noise of the 2-d spectrum 
we are not able to test whether our measurements of the broad 
Ha feature could be corrupted by a possible rotating disk as 
in the case of M84 (e.g., Walsh, Barth, & Sarzi 2010); the 
match in the widths of the [S II] and [N II] narrow features 
(as long as that of the narrow Ha) argue however against such 
a significant effect. 



,/blend 
ftia 

FWHM(Ha'""'""' 
Av 

log FiHa'""'"^) 
log L(Ha''"'"i) 



0.80 
0.93 
1570 
140 
-13.38 
37.96 



fraction of broad Ha to Ha4-[N II] blend, 
fraction of broad Ha to total Ha 
± 110; inkms-' 

± 20; in km s"' ; broad relative to narrow Ha 
observed flux (erg s"' cm"~) 
observed luminosity (erg s"')^ 



" assuming distance d = 4.3 Mpc. 

Measurements related to the broad Ha feature from the 
STIS measurements are listed in Table [T] The frac- 
tional contribution of this broad component to the to- 
tal flux of the HaH-[N II] blend is 80%. Compared to 
other galaxies with broad line emission in the Palomar sur- 
vey dHo. Filippenko & Sargent! Il997bl) . the width of the 
broad Ha is typical (1570 km s~^). However, at just 
9.1x10"'^ ergs s"', the Ha luminosity is lower than any 
known broad line source s in the Palomar sample or the 
Sloan Digital S ky Survey dHo. Fihppenko & Sargent 111 997bl: 
iGreene & Ho Il2007i) . This includes the famous low luminos- 
ity Seyfert 1 galaxy , NGC 4395, which has a lu minosity of 
1.2 X 10^^ ergs s"' (Filippen ko & Sargent I I 1989^) and is at a 
similar distance to NGC 4736. Note also that the broad Ha 
feature is redshifted by 140 ± 20 km s"' relative to the nar- 
row emission, that is locked at the systemic velocity of the ob- 
ject; velocity of fsets between the broad Ha and n arrow lines 
are not unusual ( Ho. Filippenko & Sargent II 1997al) . 

Evidence for broad Ha emission has also been presented 
in a recent PCA to mography study applied to this nucleus 
(ISteiner et al. II2009I) . The observations reported in this case 
come from the Gemini Multi Object Spectropraph (GMOS)- 
IFU data cube, and have been obtained 4 years after the HST- 
STIS spectrum had been acquired. Thus, the broad Ha line 
associated with the nucleus of M94 appears to be persistent 
for at least 4 years. They derive a broad Ha luminosity of 
6 X lO-*^ ergs s"', > 6 times brighter in the GMOS data than 
the one we detect in the STIS spectrum. The difference could 
result from true variability in the broad line luminosity, or 
could be due to d ifferences in aperture size or placement. We 
note that lSteine r et al. (2009) suggest the broad line region is 
offset from the photocenter of the galaxy by 0"15, while our 
spectral 0."25 x 0"! aperture coincides with the photocenter 
Even with their higher Ha luminosity, this source is among 
the least luminous known broad line AGN. 

2.2. Ground-based optical spectroscopy & the spectral 
classification 

Ground-based optical spectra of the NGC 4736 nucleus 
are available fr om the Palomar spectroscopic survey by 
IHo. Filippen ko & Sargent I (|1997al) and the integrated spec- 
trophotometric survey survey of iMoustakas & Kennicutt I 
(i2006.) conducted with the 2.3m Bok telescope, as well as 
from new data we acquired on February 2008 at the MMT 
Observatory. The MMT spectrum is obtained with the Blue 
Chanel Spectrograph, the 500 grooves/mm grating used in 
first order with the l"slit, and covers AA3800-7000 with 
3.6A resolution. None of these data show indications of broad 
Ha emission. This outcome is not surprising given the fea- 
ture's flux in the spectrum, which would be very difficult 
to discern in the l"or 2"apertures employed in these ground- 
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Fig. 2. — NGC 4736 within emission-line diagnostic diagrams. Filled cir- 
cles reflect measurements b ased on the Palomar obse rvations. Diamonds 
represent measurements fromlMoustakas & Kennicutt I E006), while crosses 
show measurements from our MMT spectrum. The vertical bands indicate 
the HST-STIS measurements and the associated errors (with no coirection 
for reddening; [O III] and H/3 are not available in the HST spectra). The solid 
(black) curves indicate the Kewley et al. (2006) classification, while the dot- 
ted (blue) lines indicate criteria used bv lHo . Filippenko & Sar gent i tl997aft . 
The backgro und grey points correspond to m easurements of SDSS nearby 
galaxies from lConstantin. Hovle. & Vogelevl>2008i) . 



TABLE 2 
Emission Line Measurements' 



line name 


HST-STIS 


Palomar*" 


Bok"' 


MMTO'' 


[0 I] A6300 


5.03±0.31 


6.03 


2.8±0.9 


7.0±1.0 


Ho (naiTow) 


3.00±0.32 


> 25.12 


13.5±1.9 


22.0±2.4 


[N 11] A6583 


5.25±0.25 


54.01 


21.5±1.5 


39.3±2.9 


[S 11] A6716 


2.48±0.18 


18.59 


8.3±1.1 


15.5±0.3 


[S 11] A6731 


2.16±0.18 


16.33 


7.2±1.1 


9.8±0.3 



" All fluxes are in units of 10 erg s ' cm ^, and represent the 
observed values, not corrected for reddening, 
nonphotometric conditions; line ratios are at least 50% uncertain. 
2"5 slit, 2.3m teles cope; photometric conditions 
dMoustakas & Kennicutt 120061) . 
" l"slit; fluxes are in units of 10"" erg s"' cm"^. 



based observations, which are at least one order of magnitude 
larger than that used in the HST-STIS observations. Through 
a 2"aperture 40 pc), NGC 4736's emission complex Ha+ 
[N II] is generally heavily swamped by the host stellar light. 

Probably as expected, measurements of NGC 4736's 
narrow-line emission are consistent with this system being 
powered at least partially by a non-stellar source, regard- 
less of the resolution of its observation. Figure |2] shows 
the location of this nucleus in the 3-dimensional diag- 
nostic diagram usually employed in classifying em i ssion- 
line galaxi es (e.g., Baldwin, Phillips, & Terlevich' 119811; 
IVeilleux & O sterbrock 1987; Kewley et al 2006i), for 
different sets of data mapping ~ 10 to ~ 100 pc. As with 
the HST-STIS observations, the ground-based measure- 
ments are performed after the host stellar contribution has 
been subtracted from the observed spectrum by means of 
absorption galaxy template fits. Our MMT spectrum offers 
only an upper limit for the [O III]/H/3 while the Palomar 
line flux measurements are at least 50% uncertain. We list 
the measurements of the nebular emission of this nucleus as 
observed by Palomar, Bok, MMTO, and HST-STIS in Table 
|2l Because measurements of the [O III]/H/? ratio are not 
available in the HST spectrum, we show the high-resolution 
measurements in Figure |2] using the [O III]/H/3 from the 
Palomar catalog; it is readily apparent that the high resolution 
data are consistent with a LINER or a Seyfert classification 
for this object, and reveal its type 1 (broad-lined) character 



2.3. The radio, UV and X-ray observations & the astrometric 
coincidence with the broad Ha emission 

Across the electromagnetic spectrum, the center of this 
galaxy exhibits a complex morphology. There is a plethora of 
bright radio, UV and X-ray sources detected in the center of 
NGC 4736. There is a nuc lear compact (15 GHz) radio source 
measured bv lNagar et al. 1 (2002) that appears to be associated 
with that of the brightest of the two close (8.49 GHz) radio 
sources revealed by Kording et al. (2005J, and to that of the 
brightest of the two u ltraviolet point sou rces, which varies on 
a 10- year timescale dMaozetal. 112005b . Chandra observa- 
tions dEracleous et al. II2002I) reveal numerous discrete X-ray 
sources in the inner galaxy; the second brightest X-ray source, 
X2 (Lx.2-mev = 5.9 x lO"'^ erg s~'), coincides within the er- 
rors with the nucleus position. The off-nuclear radio, UV and 
X-ray sources are apparently unrelated to each other 

The relation to possible optical counterparts of these ob- 
servations is not well constrained in the literature. We thus 
investigated this issue, with a particular interest in the degree 
to which our newly detected broad Ha coincides in position 
with the multi-A detections. We used for this purpose a vari- 
ety of archival images from the HST: WFPC2/PC data in the 
F555W filter data (PID: 5741 and 10402), HRC data in the 
F250W and F330W filters (PID: 9454), and NICM0S/NIC3 
data in the F160W filter (PID: 9360). Afl data were down- 
loaded from the HST archive and images were drizzled to- 
gether when required. Absolute astrometry was performed on 
these data as well as the STIS observations (taken with the 
slit out) using a ground based V-band image obtained from 
the Spitzer Infrared Nearby Galaxies Survey (SINGS) ancil- 
lary data (jKennicutt et al. 2003). These SINGS observations 
were aligned to the USNO-B system using ^80 stars and then 
astrometry of the HST data was obtained by degrading the res- 
olution of the F555W image to match the SINGS image. All 
other //5r images were then matched with the astrometry cor- 
rected F555W frame to an accuracy of < 0"05. The absolute 
error on this astrometry is about 0."2, and is dominated by 
scatter of stars in the SINGS image relative to the USNO-B 
positions. 

Figure [3] illustrates the result of this data compilation and 
the corresponding radio, UV and X-ray source matches. It 
is clear that, within < 0."2 (4 pc), there is an obvious astro- 
metric match in the nuclear X-ray, UV, optical, radio compact 
sources, and the newly detected broad Ha emission line. The 
nucleus position as observed by STIS has RA: 12''50'"53?20, 
DEC: 41° 07' 13 "40. The off-nuclear XI, that is the bright- 
est compact X-ray source detected in this galaxy nucleus, has 
no counterpart at other wavelengths. Same is true for the 
off-nuclear 8.49 GHz detection, that is only l"(^20pc) away 
from the nucleus. The off-nuclear UV source is the only one 
detected in optical light. 

3. NGC 4736 S NUCLEAR EMISSION ACROSS THE 
ELECTROMAGNETIC SPECTRUM AND Lbol 

With the wealth of data available for this galaxy nucleus, we 
are able to build its least contaminated X-ray to radio nuclear 
spectral energy distribution (SED). The multiwavelength ob- 
servations of the sources detected in the very central regions 
of NGC 4736, plotted in uLi' units, are displayed in Figure H] 
The X-ray detection X2 is represented here by a power-law, 
estimated based on its photon-indices F = 1.6, where the ab- 
sorbing column density is fixed at the Galactic value. Both 
F250W and F330W UV observations are plotted, only for the 
nuclear detections. The optical data are represented by the 
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Fig. 3. — The mosaic-ed F555W (V-band; combination of 5 WFPC2 images taken at two different epochs) image of the NGC4736 nucleus, with intensity 
contours overlaid (green). The position of the STIS-nucleus is indicated (by the yellow slit) along with the locations of the 2 most l uminous and closest to the 
nucleus hard X-ray sources, XI and X2, as detected with Chandra (red squares), along with the position of the 15 GHz VLA radio core i Nagar et al.l2 002. orange 
circle), and those of the two radio compact sources detected by Kording et al. 1 2005) (white circles). The blue and cyan contours con'espond to the F330W and 
F250W observations from lMaoz et al. 1 12005) ; note the presence of NGC 4736b, ^ 2" 5 to the north of the nucleus. The absolute astrometry is good to within 
0;'2 (pixel size is 0^05 for F555W, and 0;'025 for F250W and F330W). 
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Fig. 4.— Radio to X-ray n uclear SED of NGC 4736 superposed on SEDs 
of LLAGN from IHol 119991) . The optical STIS spectrum stands out as an 
amorphous blob featuring the strong Ha emission feature. The 2-10 keV 
Chandra X-ray detection X2 is depicted as a power-law corresponding to 
r = 1.6. The low spatial resolution Spitzer IRS observations are shown for 
comparison. 

HST -STIS spectrum, featuring the strong Ha emission Une. 
Measurements corresponding to all the radio observations dis- 
cussed above are indicated. For the sake of completeness, we 
also include in th is plot lower spatial re solution observations 
from SINGS IRS jKennicutt et al. 12003b as well. We show all 
of th ese m easurements superposed on the SEDs of LLAGN 
from iHo 1 (01999). No artificial normalization has been per- 
formed, and there is no correction for absorption in either the 
Galactic extinction (except for the X-ray data), or intrinsic to 
NGC 4736. 

With X2 and n4736-b as the X-ray and the 8.5 GHz coun- 
terparts of the broad Ha detection respectively, and assum- 
ing the continuum could be described as simple power-laws 
between the points we present data for, this object's nuclear 



TABLE 3 
Nuclear SED Data 



log (u/Uz) 


log {vLvlergh) 


Resolution 


Obs. date 


Instr. 


18.38 


38.69 


0"15 


5/2000 


Chandra-ACIS 


17.68 


38.42 


0"15 


5/2000 


Chandra-ACIS 


15.08 


39.41 


0"5 


6/2003 


ACS F250W 


14.96 


39.72 


0"5 


6/2003 


ACS F330W 


14.68 


39.91 


0"1 


6/2002 


HST-STIS 


14.64 


39.93 


0"1 


6/2002 


HST-STIS 


13.48 


40.58 


4;'7 


5/2004 


IRS-SH 


13.19 


40.33 


11"1 


5/2004 


IRS-SH/LH 


12.91 


40.87 


11"1 


5/2004 


IRS-LH 


10.18 


35.79 


0"15 


1/2001 


VLA, A config. 


9.93 


35.47 


0"24 


6-10/2003 


VLA, A config. 



SED corresponds to Lboi ~ 2.5 x lO'*" erg s"'. Spitzer data 
is not included in estimating Lboi because they do not reflect 
the IR emission of the uncontaminated nucleus; the aperture 
used in these observations is very large (see Table [3]), includ- 
ing the entire field-of-view shown in Figure [3] This luminos- 
ity makes NGC 4736 one of the least luminous LINERs with 
strong evidence of BH accretion. This source is thus a crit- 
ical signpost of BH accretion at extremely low levels. Note 
that this obj ect's weak emission is most likely not caused by 
obscuration: iHo. Filippenko & Sargent] (Il997al) list a Balmer 
decrement of Ha/ HP = 3.1 (al beit highly uncertain, w ith a 
probable error of ±100%), and lEracleous et alTI (120021) pro- 
vide an upper limit for the neutral hydrogen column density 
of Nh < 3.3 X 102"cm-2. 

3.1. Comparison with other observed LLAGN 

Figure 4 shows that NGC 4 736 's SED is very similar to 
those of other LLAGN dHo 11999.) . This suggests that, in spite 
of the significant difference in Lboh there is no fundamental 
transition in the accretion mode in this source compared to 
other LLAGN. 
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The similarity of the SED means that previously proposed 
bolometric corrections for LLAGN appear to work quite well 
for NGC 4736. The value estim ated based on th e correction 
to the [O III] line luminosity that lHeckman et ai~l (2004) pro- 
posed to work well for Seyferts is L[,oi ~ 3500 x L[oiii] = 4.5 x 
10'*° erg s"' ( here L[oiii] is measured from the ground within 
the Palomar survey, and is not corrected for intrinsic redden- 
ing). If the more recent assessment of the Lb oi = 600 x Lro im 
bolometric correction of Kauffmann & Hec kman"! (120091) is 
used, with L[oiii] corrected for reddening based on the Balmer 
decrement listed above and a r oc A""'^ attenuation law 
dCharlot & Fall 2000), then Lboi = 0.93 x 10"*" erg s"', which 
compares satisfactorily with the measured Lboi- The average 
bolometric correcti on to t he o bserved 2-10 keV X-ray lumi- 
nosity proposed by iHo I ("2008), with L/Lx = 16, results in 
a somewhat lower value, 0.94 x lO"*" erg s~'; with the more 
recent L/Lx = 50 av erage bolometric correction to the 2 - 
10 keV luminosity of lEracleous. Hwang. & Flohic I (l2010al) . 
Lboi = 1 .8 x 10"*" erg s"', and thus, very close to our integrated 
value. Within a typical uncertainty of '--^50%, all of these es- 
timates are consistent with the measured Lboi value. Other, 
more uncertain bolometric indicators are also consistent with 
our integration; the correlation between the mid-IR [Ne V] 
e mission-line lum inosity and bolometric luminosity derived 
bv lSatvapal et al. I (12007) for a small sample of much brighter 
nearby AGN, that has a large scatter 1 order of magni- 
tude), gives a Lboi = 1.6 x 10'*'^ erg s"', which is very close to 
our measured value. 

This good match among quite a variety of bolometric esti- 
mators derived independently from the multiwavelength prop- 
erties of AGN supports the AGN interpretation for this sys- 
tem. However, there are two details of its emission spectrum 
which are unlike typical broad-lined LLAGN. 

First, despite that fact that the broad Ha emission of 
NGC 4736 is one of the weakest measured among type 1 
AGN, the total Ha emission is peculiarly strong relative to 
the X-ray counterpart. The ratio Lx(2-10 keV)/LHQ is '--^6, 
and thus lower than the median of ^ 15 ex hibited by the 
type 1 AGN (and low-z quasars) included in le.g.. Ho et al. I 
(120011) study that revealed a relatively tight correlation be- 
tween the two types of emissi on. With the higher broad Ha 
flux from lSteiner et al.1 (|2009), the ratio is even lower, Lxi2- 
10 keV)/LHa 1. This finding is surprising given that the low 
luminosity nearby galaxies that deviate from the Lx-Lua lin- 
ear scaling toward lower values of the Lx/Lua ratios are the 
type 2 sources, mostly Transition Objects, where the ioniza- 
tion mechanism is not necessarily dominated by an AGN-type 
of source. Following e.g.. Ho et al. (2001) arguments, the un- 
usually low Lx/Liia ratio measured in the nucleus of NGC 
4736 suggests that either: 1) the optical line emission is not 
powered exclusively by a central AGN, or 2) the X-ray emis- 
sion in this system arises, at least partially, from a non-AGN 
source, e.g., an X-ray binary. A final possibility is that the 
source is highly variable (as suggested by the UV, X-ray, and 
Ha observations) and thus, the unusual Lx/Lua is the result 
of measurements made at different times. 

Second, NGC 4736's nebular emission shows a num- 
ber of peculiarities. To start with, the electron density in 
NGC 4736's nebular emission is very low, and it does not 
show the typical gradient exhibited by AGN, i.e., increase 
toward the more nuclear regions (Constantin et al. 2011, in 
prep.). In both the Palomar and HST -STIS observations, the 
ratio [S II] A A 6716/6731 is - 1.26 (Table |2]i, implying an 



electron density of tie ^ 10^ cm ^, which lies at the lowest 
end in the distri bution of particl e densities measured in the 
Palomar objects dHo et al. ll2003l) . 

Given the lack of a density gradient in the line emitting 
region it is then rather peculiar to observe a significant in- 
crease (by a factor of ~6) in the [O I]/Ha line flux ratio in 
the HST spectrum relative to the large aperture observations, 
which would usually be interpreted as an indication for a more 
pronounced AGN-like ionization in the more central regions. 
There are two reasons for this: 1) [01] requires a significantly 
hard radiation field, i.e., that of an AGN, to sustain a suf- 
ficiently extensive partially ionized zone in clouds optically 
thick to Lyman continuum, and thus to produce a strong such 
feature. Since the ionization potential of [O I] matches that of 
H very well, large differences in the [O I] /Ha ratios are ex- 
pected between accretion and non-accretion sources. 2) given 
the lack of a density stratification in this nucleus, the degree of 
ionization of the emitting gas is expected to diminish with ra- 
dius as the ionizing radiation emerging from a nuclear source 
falls off in density as r"^, and thus produce strong gradients in 
the [O I] /Ha ratio. 

In summary, NGC 4736 may be a broad-lined AGN, with an 
SED similar to other LLAGN, but it is atypical in several as- 
pects whose physical origin remains unclear. Given the above 
listed possible explanations for the peculiar Lx /Lua ratio, we 
evaluate in the next subsection whether photoionization by the 
weak AGN in this system is sufficiently powerful to balance 
the emission cooling in this system. 

3.2. Comparison of ionizing and emission-line power 

The multi-wavelength observations of this system allow for 
a relatively rigorous assessment of whether the photoioniza- 
tion by this system's AGN can power the measured emission- 
line luminosities, a nd in particular that of the bro ad Ha com- 
ponent. Following lEracleous. Hwang. & Flohic I (I20i0b). we 
can run an energy budget test via a direct comparison of 
the Ha luminosity and count rate with the ionizing lumi- 
nosity Li = LiRy—iookeV and the ionizing photon rate Qj = 

GlRy— lOOkeV- 

It is important to treat separately the broad and narrow 
emission line features as there is strong evidence that they 
originate from regions of significantly different physical con- 
ditions. The broad Ha comes from a much more compact and 
much denser emitting gas than the narrow Balmer and for- 
bidden lines. The difference in density is at least 3 orders of 
magnitude; the critical density for collisional excitation of the 
[O I] (~ 2 X 10^ cm"-') can be used to estimate the gas den- 
sity of the broad line region, as [O I] does not exhibit a broad 
component in this nucleus' emission. The difference in size 
is expected to be ~ 2 orders of magnitude (e.g., reverberation 
mapping of nearby Seyfert galaxies measured to be < 1 week, 
Denney et al. 2010). Thus, the mechanisms that can operate 
in these two regions are expected to be qualitatively differ- 
ent. The photon and energy balance conditions should reflect 
these differences, and thus should differ as well; in particular, 
potential contributions to the Ha emission via collisional ex- 
citation should be minimal for the narrow line emitting region, 
but important for the broad component. 

Based on photoionization models (i.e.. Cloudy, v94.0; see 
Feriand et al. 1998) computed by iLewis et al.1 (l2003h for 
a wide range of ionization parameters, densi ties and metal- 
licities, lEracleous. Hwang. & Flohicl ( 1201 Obi) finds that en- 
ergy balance in a line-emitting nebula requires that L, > 
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TABLE 4 

Ionizing and Emission-Line Power 





log Llia 


' log Ghcv" 






Palomar 


> 37.7 


>49.3 


<93 


< 2.75 


Bok 


37.5 


49.0 


165.9 


4.9 


MMTO 


37.7 


49.2 


104.7 


3.1 


HST-STIS, narrow 


36.8 


48.4 


759 


22.4 


HST-STIS, broad 


37.9 


49.5 


57.5 


1.6 


PCA tomography 


38.8 


50.3 


8.7 


0.25 


^ the luminosities are 


measured 


in erg s ' and the photon rates in s ' . 



the minimum energy balance condition for AGN ionization is given 
by Lj/Luc 180 (when a fraction of 10% of ionizing photons photons 
is absorbed by the line-emitting gas) 

' the minimum photon balance condition for AGN ionization is given 
by Qi/QHa > 2.2 for the narrow line regions, and Qi/Qua 0.25 for 
the broad component. 

18(±2)LHc(/./e, where fc is the covering factor, or the fraction 
of the ionizing luminosity of the AGN that is absorbed by the 
Hne-emitting gas. These models are covering electron densi- 
ties that appear to encompass both the broad and the narrow 
line emitting regions in NGC 4736; thus, with only a frac- 
tion of 10% of ionizing photons photons being absorbed by 
the line-emitting gas, the minimum energy balance condition 
for AGN ionization is given by L,/Lhq > 180, for both nar- 
row and broad emission features. In the same time however, a 
minimum requirement for photon balance can be quite differ- 
ent for the two emission regions: for the narrow line compo- 
nent Qi > 2.2 Qho, corresponding to the case B recombina- 
tion (i.e., one Ha photon is emitted for every 2.2 recombina- 
tions); for the denser broad line emitting region, the number 
of H-alpha photons that can be produced for each ionization 
can be at least 7-8 times high er than the standa rd case B esti- 
mate, or Qi > 0.25 Oh^ (e.g.. lOsterbrock llT989i) . 

With the Lj and Qj v alues already calculated by 
lEracleous. Hwang. & Flohic I (2010b) by integrating M94's 
nuclear SED assuming that pairs of points could be connected 
by a power-law, and with the Ha measurements from all HST 
and ground-based observations presented above, we can pro- 
ceed with the comparison. Table |4] lists the Lhq, Qua, and 
associated ratios Lj/Lna, Qi/ Qua for all of these optical spec- 
troscopic observations, where the narrow and broad Ha mea- 
surements are shown separately. 

It is readily apparent that while Li/Lua is in general > 18 
(fc = 1), it is almost never > 180 (fc = 0.1); the only excep- 
tion is the narrow Ha emission measured in the HST aper- 
ture. Thus, a dominant AGN ionization of the narrow line 
region in this nucleus is definitely possible in the HST aper- 
ture, but only for fc > 18% at larger radial distances. For the 
HST broad Ha feature, an AGN ionization is possible only if 
fc > 30%. The PCA tomography measurement of the broad 
Ha argues, however, against a balanced energy budget origi- 
nating entirely in an AGN-like power mechanism, even when 
a maximum covering factor is considered; the AGN-produced 
energy falls short of the required amount by at least 50%. In- 
terestingly, the measured Qi/Qua ratio is well within the re- 
quired photon balance corresponding to an AGN excitation 
for both the narrow and the broad emission features. Thus, 
with one clear exception, the AGN in NGC 4736 appears to 
be capable of providing enough photons to explain the ob- 
served Ha luminosity but only for relatively high f - values. 

Simply because the AGN-like SED of NGC 4736 can ex- 
plain the majority of the ionization energy and photon rate 
does not imply that the actual mechanism is an AGN; alter- 



native excitation mechanisms must be explored. Possible op- 
tions are: i) we are missing ionizing photons from accretion 
onto the central BH, or ii) there are other power sources that 
could make up the power deficit in this system, particularly 
for producing the broad Ha feature. 

The first alternative could be possible if we were observing 
an "echo" of a previous epoch of more violent accretion, a 
few hundred years ago. This idea has been explored by Era- 
cleous et al. (1995), who showed that the reverberation of an 
ionizing flare in the nebula can produce LINER-like emission- 
line ratios. In this scenario, it is expected that the central UV 
source and the [O III] line would also follow the decay of the 
ionizing continuum; while the UV observations are not pro- 
viding clear evidence for such a decay over the course of one 
decade, the multiple optical spectroscopic observations that 
we present in Figure |2]are consistent with a possible decrease 
in time (a few years) in the [O III] flux (in this model, H/3, 
[S II], [O I] are expected to decay very slowly, in 60-250 
years). Nevertheless, this duty-cycle hypothesis also requires 
that the broad-line emission fades immediately if the ionizing 
continuum declines, and the PCA tomography observations 
show that this is not happening in this source. The echo of 
such an ionizing continuum flare should also be detectable 
in a narrow-band [O III] A5007 image in the form of a ring, 
which is not readily observed in M94. 

The second alternative of power sources other than accre- 
tion, has been often proposed in explaining the emission- 
line spectra of LINERs, with mixed success (see §1). The 
most probable alternative sources appear to be the mechanical 
power delivered by compact radio jets, along with photoion- 
ization by young or post-AGB stars from old or intermediate- 
age stellar populations. Shock models are highly unlikely to 
produce broad emission features, and thus, are not favored in 
this case. Recent star-formation activity remains however a 
viable options, particularly in light of relatively new discover- 
ies of peculiar supernovae with broad Ha features that do not 
appear to fade in time. This idea is discussed in more detail in 
§5. 

4. THE BLACK HOLE MASS 

Given the unusual energetics of M94's nucleus, it is impor- 
tant to investigate whether the BH mass estimators derived for 
rapidly accreting Seyferts, which appear to be widely used for 
AGN, also work in the low-luminosity regime flagged by this 
particular system. 

There are a variety of indirect methods that can be used to 
estimate the mass of BHs in galaxy centers. The available 
multi-wavelength measurements of the nuclear emission for 
NGC 4736 allow the calculation of Mbh based on four differ- 
ent techniques, as well as with the M-a* relation. In this sec- 
tion, we explore and compare the results of these five meth- 
ods, along with their consequences for this object's energet- 
ics, and a comparison with a recent dynamical measurement 
of 6.68(5.14-8.22) x lO^'M^^ for this o bject, which is listed in 
iKormendv. Bender. & Cornell I (1201 Ih as obtained from Geb- 
hardt et al. (201 1, in preparation). We first present the meth- 
ods and associated Mbh calculations, and then discuss the 
shortcomings of each measurement. 

1. Using the M-a* relation, established for quiescent 
nearby galaxies, including both ellipti cals and spirals 
with classical bulges, as quantified by iGult ekin et al J 
( I2009 ai for ct* = 110±5 km s"' dBarth et al. ii2002l) . 
Mbh = 1.05 ±0.64 x lO^M©; within errors, this value 
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agrees well with Gebhardt et al. dynamical measure- 
ment. With this value, Lboi/i'Edd ~ 2 x 10~^, which is 
consistent with the range of values within LINERs are 
expected to lie (e.g.. Ho 2004; Ho 2008). 

2. Using the HST measurements of the FWHM and the 
luminosity of the broad Ha component (see Table [U 
within the scaling relation based solely on observa- 
tions of this broad emission feature, as derived by 
lGreene&Ho] (^05). we obtain Mbh ~ 3 x lO'^M©; the 
fractional uncertainty associated with this measurement 
is '--^30%, and includes both the scatter in the scaling re- 
lation and the errors in the line measurements. With the 
broad Ha lin e luminosity m easured in the PCA tomog- 
raphy study (Stein er et al. II2009I) . the BH mass would 
mcrease (by a factor of 6.5) to Mbh = 1 .9 x lO^M© . The 
difference in these two values could be considered the 
most conservative uncertainty associated with this BH 
mass estimate. The corresponding Eddington ratios for 
these BH masses are Lboi/iEdd « 1 - 7 x lO""*. 

3. Mid IR detection and measurements of the [Ne V] 
14.32 ^m) and [O IV] ( 25.89 /im) emission lines 
Dudik. Satvapal & Marcu l l2009) give, via the empir- 
ical correlations between the MIR line luminosities and 
reverberation rnapping-based Mbh values presented by 
iDasyra et aTl j20()8h . a black hole mass of 2.3 ± 0.4 x 
IO^Mq, and 1.7 ±0.5 x IO^Mq respectively. The cor- 
responding Eddington ratio is in this case ^^boi/^Edd ~ 
1.0 X 10"^. Note that the scatter adopted for these re- 
lations is only a lower limit of the real value, thus, the 
uncertainty may be larger. 

4. The norm aUzed X-ray excess var iance method, as de- 
scribed in Papadakis et al. (2008), applied to Chandra 
observations (Eracleous et al. 2002) gives Mbh = 2.5 ± 
1.7 X IO^Mq. To be specific, we used for this calcula- 
tion the excess variance a = 0.06 ±0.04, Lboi = 2.5 x 
lO"*" erg s"' as derived in Section [3] and v\f= \/T = 
1/14 h"' (where T is the length of the light curve), 
to estimate the break frequency Vhi and then Mrh 
via eq uations 4 and 6 respectively, of iPapadakis et al.1 
(I2008IV 

5. The Mbh of this system can also be obtained via the 
"fundamental plane of black hole activity" that re- 
lates black hole mass to the emitted compact radio 
Lr = vL^iSGHz) and hard X-ray luminosities Lx(2- 
10 keV), and spans nine orders of magnitude in black 
hole mass. We calculate Lr = 1.7 x 10-^^ erg s"' us- 
ing the two 8.5GHz and 15 GHz radio measurements 
presented in this paper (Figure [3] Table [3]), with a flux 
modeled b y a power-law Si, oc v'"^ . Us i ng the empir- 
ical fits of lMerloni. Heinz & di Matteol (l2003h for the 
fundamental plane relation, log (Mbh/Mq) = 5.9 ± 1.1 
for the black hole in this system; note that this mea- 
surement and its associated errors embrace all of the 
estimates presented above, and thus, do not provide 
any additional constraint to the Mbh- Interestingly, 
the latest derivation of the fundamental plane relation 
(iGiiltekin et al. I l2009bl) . applied to very low nuclear 
galactic luminosities, provides a log (Mbh/Mq) = 7.2 ± 
0.4 for this system, which remains consistent only with 



the M — a* value, as it departs considerably from those 
provided by the relations employing AGN emission. 

There is a significant inconsistency between the value given 
by the M-cr* relation, which is supported by Gebhardt's dy- 
namical measurement, and those based on the AGN emis- 
sion. The M-a* relation suggests a BH mass in NGC 4736 
of 1x10^ M0, and a correspondingly low Lboi/i^Edd of ^ 
2 X 10"^. Three of the estimates based on the AGN emission 
converge to black hole masses of ^ 10'' M0 (and Lboi/^Edd ^ 
10"^), showing a surprisingly consistent departure of two 
orders of magnitude from the M-cr* estimate. The fun- 
damental plane relation provides a Mbh value right in be- 
tween these two different situations, however, with no real 
additional constraint, due to its associated large uncertainty. 
These differences are somewhat puzzling given that the BH 
mass estimates based on nuclear emission properties are all 
calibrated to follow the M-a* for high mass black holes 
(Mbh ^ IO^Mq). It is however true that the calibrators are 
biased toward nearby Seyfert galaxies with much higher Ed- 
dington ratios than that of NGC 4736. We briefly discuss in 
the following subsections more specific weaknesses of each 
measurement. 

4.1. Caveats of the M — a* relation 

The M-a* relation is expected to provide a reliable es- 
timate of the BH mass as it is based on the strong corre- 
lation between dynamical mass measurements of supermas- 
sive BHs and their hos t properties j Ferrarese & Meritt 2003 
Gebhardt etani2000t iTremaine et al. 2002; Gulteki n et all 
'2009a|). This relationship is derived primarily from elUp- 
ticals and spirals with classical bulges (formed during ma- 
jor mergers). Recent observations suggest it may not be 
valid for samples of later-type spirals which more commonly 
host pseudobulges (formed via secular disk processes) (e.g. 
iGreene et al. II2OI0I) . Because the distinction between classi- 
cal and pseudobulges is based on for mation, it does not simply 
corre l ate with observable p roperties ( iKormendv & Kennicutt I 
120041; iFisher & Drorv 1 120 10). In general pseudobulges are 
less luminous, have lower bulge-to-total ratios, have on-going 
star formation and lower Sersic indices than classical bulges. 
The possibility that the NGC 4736 is in fact a pseudobulge 
provides a solution to the apparent conflict in BH mass esti- 
mates. 

Nevertheless, pseudobulges are difficult to identify, and 
there is not yet a consensus as to what defines them. To 
complicate things further, classical and pseudobulges can ex- 
ist within the same galaxy (e.g., NGC 2787; Erwin et al. 
2003). The classification of th e NGC 4736 bulge is am- 
biguous: iFisher & Drorv I (120101) classifies it as a pseudob- 
ulge based on its nuclear spiral and bar and low Sersic index 
(n = 1.3). However, they also find it has a low star formation 
rate, more typical of classical bulges, and thus classify it as an 
"Inactive Pseudobulge;" if NGC4736 hosts a pseudobulge, it 
is an atypical one. We have created a surface brightness pro- 
file from NICMOS and 2MASS H-band data from the Large 
Galaxy Atlas (Jarrett et a l.ni2003h and fou nd results that con- 
flict with the fits of Fisher & Drorv I (1201 0). Specifically, we 
find Sersic indices for the bulge of « = 2.3-3.0 depending 
on the radial range and type of fit (single vs. double Ser- 
sic), which are consis tent with a classical one (as shown in 
IFisher & DrorvllIoToh . 

Our surface brightness profile fits also reveal the pres- 
ence of a nuclear star cluster within the central ~ 0.6" 
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(~12pc), with an H-band magnitude of ^12.5. Such nu- 
clea r star clusters are c ommon in early-type spiral galax- 
ies dCaroUo et al. Il2002l) . The luminosity and mass of nu- 
clear star clust ers are known to sca l e with bulge luminos- 
ity and mass (Bal cells et alTI 120031: [Ferrarese et_an 120061: 
iRossa et al. .2006.) . and the NGC 4736 nuclear star cluster has 
a lumino sity that is 0.1% o f its bulge, typical for nuclear star 
clusters jCote et al. 112006) . Nuclear star clusters commonly 
co-exist with black holes, but there are a very limi ted number 
of cases where masses for both can be estimated (iSeth et al. I 
|2008). In these cases, including the Milky Way, the BH mass 
is similar to the mass of the nuclear star cluster within an or- 
der of magnitude. There is also some evidence that the ratio of 
BH m ass to nuclear star cluster mass increases with spheroid 
mass ( Graham & Spitler II2009I) . For NGC 4736, assuming an 
old population with an H-band M/L ~ 0.7, the nuclear clus- 
ter would have a mass of ^ 2 x 10^ Mq. This mass is quite 
similar to the M-a* BH mass estimate, suggesting thus the 
presence of a similar sized BH. 

4.2. Caveats of the BLR scaling relation 

The BLR scaling relations were derived using high lumi- 
nosity systems (i.e., Seyferts, not LINERs) and have been 
scaled to match the dynamical black hole detections of 
BHs with masses > 10* Mq. Thus, their applicability to 
NGC 4736, where they yield an estimate of 10^ Mq, 
is a (perhaps unwise) extrapolation. These relations as- 
sume that the BLR is virialized due to proximity to the BH 
jPeterson & Wandel I IT999I) and thus the 5 100 A continuum 
luminosity correlates with the emissivity weighted radius of 
the BLR, and thus with the BH mass (Kaspietal. 2000). 
The overall errors associated with BH mass measurements 
based on these relations do not exceed 0.5 de x ( Vesterg aard . 
I200I iNelson et al. IllOOl lOnken et al. II200I . and there ap- 
pears to be good consistency with BH masses obtained via the 
M-a* relation for relatively bright AGN (L > lO'^^ erg s"'), 
with Mbh > 1O''M0 dBartii et al. I i2005). More recent caU- 
brations of the radius-luminosity relationship on which these 
techniques are based infer that BH masses have been overes- 
timated, however, only by up to a factor of ^3 (< Bentz et al. I 
12009). This latter study also indicates a trend toward larger 
uncertainties, and larger amount of overestimation in the BH 
with decreasing luminosity, however, no conclusive results 
are yet available for this regime. 

4.3. Caveats of the MIR line-correlation 

The iDasvra et al. I (1200 8') empirical relation between the 
MIR line emission properties and the BH mass is derived 
using reverberation mapping BH masses, and thus, as with 
the previous method, may not apply to systems with Mbh < 
1O*M0. The relation holds for systems with Lboi/^'Edd > 
0.003, but not necessarily beyond this range. The Mbh ~ 
lO^M© given by this method places this object at the low end 
of this Lboi/i'Edd range: on the other hand, a more massive and 
thus quiescent BH would correspond to a [Ne V] luminosity of 
a few orders of magnitude lower than that measured. The scat- 
ter associated with this relation remains in average 0.5 dex, 
and thus cannot account for the 2 - 3 orders of magnitude 
difference between this value and that obtained via M- cr*. 

4.4. Issues with estimating the BH mass from X-rays 

The BH mass derivation based on X-ray variability, or more 
precisely on the relation between the excess variance and 



Mbh, relies on the hypothesis of a universal power spectral 
density function (PSD) shape and amplitude in AGN, which 
is based on the idea that the X-ray variability mechanism and 
the accretion efficiency are the same for all AGN, at all red- 
shifts. These assumptions appear to hold for the obje cts in- 
volve d in deriving and testing this relation (Papadakis e t al. I 
I2OO8I) : that sample is, however, small and rather iDiased toward 
luminous (X-ray) sources, with log Lx /(ergs"') > 41 .5. None 
of the 2-10 keV sources detected in the nucleus of NGC 
4736 are in this luminosity range. Nevertheless, Galactic BHs 
(Mbh < IQ^Mq) in their hard states show variability properties 
that match well those of AGN, both of these types of sources 
falling on the same projection of the Ig - Mbh - iboi /i^Edd 
plane, where 7g = 1 /h%f, suggesting that an extrapolation to 
the intermediate mass BH (or lower L) level is practicable. 

4.5. Problems with the Fundamental Plane 

Finally, it appears that the main problem with the funda- 
mental plane is that the relation is not sharpened enough to 
provide strong constraints on BH masses. Because the rela- 
tion spans nine orders of magnitude in BH mass, it is expected 
to equally apply to any value of Mbh in the range we are in- 
terested in. It is also the case that a wide variety of BH accre- 
tion models (e.g., with efficient and inefficient flows for the 
X-ray emission, or associating X-ray flux with synchrotron 
emission near the base of a jet) are consistent with this rela- 
tion (e.g., Falcke & Biermann 1995: Heinz & Sunyaev 2003), 
suggesting that a large diversity of accretion modes or rates 
are accommodated. 

Nevertheless, the fundamental plane remains best con- 
strained only for systems with Mbh ^ IQ^Mq, and for very 
low nuclear luminosities, i.e., with negligible or zero AGN 
contri bution that allow a dyn amical measurement of their BH 
mass dGiiltekin et al. I2(j()9bh . In this regime, the latest deriva- 
tion of the fundamental plane relation provides for M94 a 
Mbh value consistent with the M-a* estimate and the dy- 
namical measurement, while it would not reliably constrain a 
^BH IO^Mq value. The scatter in the fundamental plane re- 
lation increases for lower BH mass es, with higher Eddington 
ratios dKording et ani2006; Gultekin et al. I|2009b,) . 

5. DISCUSSION: ALTERNATIVE POWER GENERATION 
MECHANISMS 

The detection of broad Ha emission, combined with the 
spatial coincidence of this emission with the detection of X- 
ray, UV, and radio compact sources, provides strong evidence 
that this galaxy hosts an accreting massive black hole. The 
SED of this nucleus makes this object one of the lowest lu- 
minosity LINER with a distinct contribution to the total emis- 
sion by black hole accretion. In this scenario, the nucleus 
is an AGN and the presence of off-nuclear sources, particu- 
larly the radio and the UV detections, may result from rem- 
nant jet activity emerging from the nucleus. It is definitely 
exciting to detect AGN activity at energy levels equal to that 
of several young su pernova remnants of the Cas A variety 
dTurner & Holll994|). an O B associ ation which hosts a high 
mass X-ray binary dEracle ous et al . 1120021), o r simply a group 
of five late O supergiants (Ma oz et al. 111995b . 

Nevertheless, the presence of a number of unusual off- 
nuclear sources, coupled with the apparent deficit in the pho- 
toionizing photons, and with the fact that the BH mass es- 
timates based on AGN emission appear to fail for this ob- 
ject, encourage us to explore alternative scenarios for the 
NGC 4736 nucleus. In particular, it is very likely that the BH 
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mass estimates that exploit the multi-wavelength AGN char- 
acteristics don't work for this system because at least some of 
this emission is not the result of BH accretion. 

There are certain (peculiar) kinds of core-collapse SNe 
that present multi-wavelength observations, and in particular 
broad Ha emission components with characteristics that are 
very similar to those we measure in the nucle us of M94. The 
SN 2005ip, presented bv Smith et al. I (120091) . is a very good 
example. The so called intermediate Ha component associ- 
ated with this SN ejecta presents the same FWHM and bright- 
ness level as the broad Ha detected in the STIS aperture, and 
it does not show any sign of diminishing its strength over more 
than 3 year period (see their Figure 7; a very broad component 
with FWHM > 10000 km s"' Uke the one exhibited by SN 
2005ip would not be measurable in the galaxy spectrum as it 
would be completely swamped in the continuum stellar light, 
even in the HST-STIS observations). 

SNe as luminous as X2, with intrinsically hard X-ray spec- 
tra (photon index T < 1), have certainly been encountered, 
for cases observed few years after the explosion (e.g., ATe 
#1023; Pooley et al. 2007); for these cases, LxlL\,o\ ratios 
are high relative to more standard SN cases (Immle r et al. I 
|2007). Another particular example of a SN that matches well 
the measurements of the nu clear emission in NG C 4736 is 
SDSS J09529.56H-214313.3 (i Komossa et aDl2()09 ). which is 
believed (but not confirmed) to be a SN type Iln. For this sys- 
tem, the broad Ha component stays strong for at least three 
years, and its Lx(2-10 keV)/LHa ~ 1, matching thus very well 
the surprisingly low value measured for NGC 4736. 

The peculiar nebular, radio and UV characteristics of the 
nucleus of NGC 4736 appear to also compare well the 
SN phenomenology. The electron densities (or [S II] line 
flux ratios) measured in this galaxy center are matching ex- 
actly the ones measured in the environments of the extraor- 
dinary type Iln SNe we compare here with, e.g., SDSS 
J09529.56H-214313.3. The two different band radio observa- 
tions of this nucleus are consistent with emission from extra- 
galactic SNe, iri both i r itensity and decline rate in flux density 
dvan Dvk et al. I [T99l l^lhams et al. I llOOl IChandra et aU 
120091) . Moreover, the compact radio off-nu clear detection, 
only 1 "(20 pc) away from the nuclear one (iKording et al. I 
120051) could be interpreted as the result of shock emission 
associated with a nuclear core-collapse SN. Measurements 
of the brightness of the off-nuclear UV detection in F250W 
and F330W place this object into the O star of late (5-ish) 
type spectral category, fitting thus well into the idea that this 
nucleus could simply be a star-forming site, and that NGC 
4736a's emission includes significant contributions from a 
SN which exploded close to the weakly active (and massive, 
Mbh lO'^-^'Mp) central BH. 

This SN contribution scenario may have its own drawbacks. 
If the broad Ha emission has actually gotten brighter by a 
factor of 6, as suggested by the comparison between the HST- 
STIS and GMOS observations, the SN interpretation becomes 



problematic; variability of the broad line region originating in 
the AGN could, in principle, account for this effect. Also, the 
observed intra-day (hour-scale) X-ray variability measured 
for X2 remains yet to be detected in a supernova, and con- 
flicts with the physical scale over which X-ray emission is 
expected in SN remnants; thus, this particular behavior may 
remain strictly associated with the AGN. Additional observa- 
tions would be necessary to fully confirm or rule out the SN 
scenario. Specifically, new high resolution observations in the 
optical, UV, or X-rays would be able to confirm if there is a 
fading in the light curve, as expected from a SN, and would 
also allow accurate localization of the source. High S/N opti- 
cal spectra would much better resolve the emission-line pro- 
file, and a possible temporal evolution. 

6. CONCLUSION 

We have presented here an exhaustive multiwavelength 
analysis of the nuclear emission properties of NGC 4736, 
prompted by new measurements of a broad Ha emission 
component detected in its high resolution HST-STIS optical 
spectrum. This broad Ha component, with a luminosity of 
9 X 10^^ ergs s"', is one of the lowest luminosity broad line 
known. This broad Ha is coincident with a compact bright 
X-ray and radio source. Our measurements of this object's 
spectral energy distribution reveal a bolometric luminosity of 
Lboi « 2.5 X 10"*" erg s"', that categorizes NGC 4736 as one 
of the least luminous LlNERs with strong evidence for BH 
accretion. Our comparison of five independent BH mass es- 
timates reveals a discrepancy of two orders of magnitude be- 
tween the value ^ lO^M© predicted by the M-a* relation 
and the value ^ IO^Mq toward which methods based on AGN 
emission activity in optical, mid-lR and X-rays, seem to con- 
verge; the fifth method is provided by the fundamental plane 
relation, which however, due to its large associated uncertain- 
ties, does not offer any additional constraint to this compari- 
son. 

We conclude that this system's BH mass cannot be reliably 
estimated via standard AGN BH mass indicators because the 
nuclear emission in this system is not entirely tracing the ac- 
cretion onto the central BH. Our assessment of the energy 
budgets of the ionizing and emission-line power suggests a 
possible deficit in the AGN ionization and production of a 
broad Ha emission feature which can be made up by a pecu- 
liar kind of Type Iln SN that matches well the nuclear emis- 
sion of NGC 4736 over the whole electromagnetic spectrum, 
and supports this galaxy nucleus' general aging starburst-like 
appearance. 

We thank the anonymous referee for constructive comments 
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work was provided by NASA through grant number HST-AR- 
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